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ABSTRACT 


This study presents the results of an experiment in objective 
analysis of oceanographic data for a limited area. The objective 
analysis is designed to provide a reliable operational system for 
tactical use in coastal waters. It is shown that this approach makes 
it possible to obtain a very detailed analysis with good vertical 
consistency and that only a relatively small amount of highly accurate 
data is required. The procedure includes a polynomial interpolation 
and a dynamical interpolation. The dynamical interpolation which is 
introduced involves the analytical solution of a finite differential 
equation by means of a rotational computing molecule. The procedure 
requires only a small computer and iittle computer time. This methed 
will provide a basis for short-time forecasts of oceanographic para- 


meters usifg only small computer centers or even time sharing systems. 
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LI. INTRODUCTION 


Analysis of oceanographic data has been traditionally a very slow 
and essentially a subjective process due to the time involved in an 
oceanographic expedition and the small amount of data obtained. In- 
creasing the number of systematic observations in time and space 
requires faster analysis procedures and objective methods. Such methods 
will provide information for regions of sparce data. The analysis pro- 
cedure should also standarize sequential analyses so that the results 
may be .compared. 

Objective analysis techniques are currently used on an operational 
basis by the Navy and other government agencies concerned with atmos- 
pheric phenomena. Background information availabie for formulating 
an objective analysis scheme for oceanographic data is primarily that 
from investigations and from experience with schemes for meteorological 
data. However, there are significant differences between both the 
type and distribution of meteorological and oceanographic data which 
have to be considered. 

The primary purpose of most objective analysis schemes is to 
transform data from observations at irregularly spaced points into data 
at points within a regularly arranged grid (grid points). A subjective 
method for doing this would be to have an analyst manually (subjectively) 
draw isolines and then interpolate (objectively/subjectively) to obtain 
values at the appropriate grid points. An example of an operational 


grid is given in Figure l. 
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FIGURE 1.- Example of the grid used by the Joint Numberical Weather 
Prediction Unit (After Cressman 1959) 








Procedures used in objective analysis have to be systematic and 
quantitative so that they can be performed by numerical computations 
on a computer. Computer processing is by its very nature, objective. 
In practice, once the observations are relayed from reporting stations 
all procedures leading to an analyzed chart should be completed within 


a minimum of manual intervention. 


A. BACKGROUND 
Research and the subsequent development of objective analysis 
methods took place in response to the need to prepare large numbers 
of meteorological observations representing the atmosphere, both in 
the horizontal and vertical, for numerical weather prediction procedures. 
Numerical weather prediction requires, in addition to values at reg- 


ularly spaced grid points, 


(D 


diting (removal of erroneous observations) 
of these data and checking for consistency and continuity in time and 
Space. These requirements would also apply to objective analysis 
schemes for oceanographic data. 

Once a grid of equally spaced values, checked for possible errors, 
is generated by an appropriate objective analysis scheme further work 
can be performed readily on a computer. These steps also could be 
performed in subjective analysis but the time required to obtain the 
completeness achieved by an objective analysis scheme would be out of 
proportion to the time needed for the computer analysis. 

An early publication on the subject of objective analysis of geo- 
physical data was by Panofsky (1949). Various approaches for per- 
forming objective analysis with meteorological data were presented 
during successive years by Gilchrist (1954), Bergthors on and Doos 


(1955), Jowmson (1957), Cressman (1959), and Smith (1902). The latter 





references describe versions of operational objective analysis schemes 
presently used in the principal weather centers of the world. However, 
most of the currently used objective analysis schemes are based on the 
"Operational objective analysis system'' described by Cressman (1959). 

The general approach in these schemes is a very straightforward 
corrective approach in which an initial (preliminary) field is modified 
by available observations and the result is then used as preliminary 
field which is again modified with the observations but with more 
resolution. The repeating corrections are called scans. This process 
may be repeated many times but in general satisfactory results are 
obtained with the first or second scan. 

In the corrective approach, an analytical model defines the cri- 
terian for checking the validity of observed values and determining 
values at grid points. Determination of grid point values trom ooser~ 
vations at Sopa stations is achieved by some form of interpolation. 
One method-of interpolating for a scalar field is to represent the 
contours of isolines surrounding the grid point as a function of 
position. This may be accomplished by means of a surface least squares 
polynomial interpolation. The order or degree of the polynomial in the 
least square interpolation depends on how detailed the final field 
should be; 

The degree of smoothing will enone upon the nature of the fore- 
cast and the particular model selected.....Smoothing to eliminate 
unwanted, though real, wave lengths has been taken to fall witha 


the province of the forecast rather than that of the analysis." 
(Johnson 1957) 


ee 


Nearby is normally considered to be within one grid distance. 


10 





Oceanographic analyses have been made only as a related factor in 
weather analysis. Climatological forecasts of ocean surface temperature 
anomalies have been made by Adem (1969), Jacob (1967), and Adem (1970) 
using the equation for the conservation of thermal energy. 

The type of data expected in an oceanographic survey introduces 
the problem of holes between transects and therefore the need of some 
type of dynamical interpolation. Typical cruise tracks appear in 
Figures 2 and 3. Provisions for lack of data for meteorological 
analysis of a given region when good data is available at the boundary 
have been studied by Thompson (1961), Richardson (1961), Smith (1961), 
and Smith (1962). The method generally used is to 'advect good data 


from outside the region (or Mhole'!)'! (Smith 1962). 


B. PURPOSE AND SCOPE OF THIS STUDY 

The purpose of this study is to develop an objective analysis 
scheme for oceanographic data which will meet operational requirements 
described ‘isa the following paragraphs. 

The high cost or operating oceanographic ships at sea is a primary 
reason to develop methods to complete instant detailed analysis. Such 
an analysis would show interesting features based on early observations 
but not considered during the initial planning which would permit 
scheduling of complementary observations during the same operation at 
sea. Furthermore, operation requirements of submarine and antisubmarine 
warfare have shown the need of faster and more accurate short term 
forecasts of oceanographic parameters which could be produced through 
the use of numerical procedures. 

The scope is to develop a scheme which would resolve fine details. 


Small thermal inhomogeneities have great effect in sound propagation 
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especially when relatively high frequencies are used. Those fine 
details cannot be observed nor forecasted with large scale schemes and 
therefore, this study considers small scale objective analysis methods 
for oceanographic data. The method will include some type of dynamical 
interpolation. The scheme will be evaluated with respect to its use 
for tactical purposes and as a "feed back" for planning during ocean- 


ographic surveys. 
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II. APPROACH 


Data obtained from the standard oceanographic stations differ from 
the meteorological data in the number of stations, the accuracy Oe 
observations and in the fact that the former nearly always has to be 
considered Eaoneice for long time steps over small areas. In general, 
an initial field will not be available in oceanographic analysis nor 
will there be enough data to obtain a complete interpolated field. 
Therefore, new procedures have to be formulated in the development of 
the objective analysis for oceanographic data. Techniques described 
above for meteorological analysis will be applicable only in a relatively 
small part in this particular problem. 


t a e o 
a numerical scheme if wnicni 2 


ids) 


The objective analysis developed i 
least squares surface polynomial is used for interpolating and smoothing 
the available data to nearby grid points. Then a dynamical inter- 
polation is used to compute values at those grid points which are in 
"holes" or regions of no data. The dynamical interpolation is achieved 


by means of a rotational molecule. 


A. POLYNOMIAL INTERPOLATION 
The first step in the scheme is to reduce the original data to 


2 
neighboring grid points using surface least squares polynomial 
i 


Synoptic data represent observations taken simultaneously. 


Zz 
Neighboring refers to data within one grid distance. 
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interpolations. A second degree polynomial is used when six or more 
data points are available in a given square around the grid point, and 
first degree polynomial when only four or five points are available. 
A weighted gradient propagation is used when only two or three data 
points are available. When only one observation is available in the 
entire square around the grid point some criteria about the distance 
is established and if this distance is, Say, less than five miles the 
best approximation will be to give the same value to the grid point, 
otherwise this value is neglected. 

The polynomial solved for the second degree polynomial interpolation 


is: 


A, + yt CL, Xd + CL, X"+ 4,4 * A, (1) 


The first degree interpolation formula is: 


ng t An X ii aX a5 Af (2) 


and the weighted eradient propagation formula iS: 


v4 4 ig W/ (z 21%, s + ¥(%- x.) (3) 
—_ (x,-x,O.-¥,) 


Where V/ = 4 0 inside the square of the grid poiat considered. 
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The matrix obtained in the first two cases is solved for Gey ; 
since the grid point is at location OO , using the subroutine SIMQ 
from the computer library of the Naval Postgraduate School. 

This scheme for obtaining values for some grid points enables the 
maximum use of data. 

In meteorological analysis the regression has been used and also 
the correction method. In the latter, it is assumed that the difference 
between the preliminary value and the observed value at a given station 
is equal to the difference between the derived value and the preliminary 
value at the grid point. Another approach in the correction method is 
to assume that the gradient of the preliminary fielavat Ene erid point 
is representative for all the area between grid points and the station 
(gradient propagation). Then the correction is a function of the 


9 


distance between grid points and station. in tnese yreceaures a 
statistically determined weight function (weighted gradient propagation) 
is required. 
In all interpolation methods, the distance between observations and 

a grid point becomes a crucial factor; for example, 

"Quite often it is not possible to get a reasonable analysis only 

by means of interpolation between synoptic observations. [It is 

quite clear that the distance between observations must be small 

compared with the size of the system to be analyzed." 

(Bergthorsson and Doos 1955) 


Therefore, distance limitations introduce the need for a dynamical 


interpolation whenever 'holes'' are present. 


B. DYNAMICAL INTERPOLATION 
The initial step in developing a suitable dynamic interpolation is 
simplification of the basic conservation equations for a scalar quantity 


in a fluid medium. In order to apply a simplified equations to a set of 


Ay 





regular horizontally and vertically distributed data, certain assump- 
tions are required about the flow characteristics. These assumptions 
have to be based on conditions observed or generally accepted for an 
oceanographic environment. The simplification of these equations, the 
assumptions involved, and the numerical procedures used to apply the 
final form of the equation are discussed in the following paragraphs. 
1. Analytical Model 

Consider the analysis of a conservative scalar field. The 
available data will represent quasi-synoptic observations provided by 
standard stations in a relatively small area of the ocean. Therefore, 
some of the following information could be available for each station: 

(1) Temperature 

(2) Salinity 

(3) Wind velocity 

(4) Dynamic height anomaly 

(5) Sigma-t 

(6) Sound velocity 
at observation depths and at standard depths. 

An analytical model for analysis of temperature, salinity and 
sound velocity would be based on the following steady state conservation 


equation for a non divergent QIEDELGLe 


GN ales 


ll Oy 20 ee -2 |) +2lka0} 2. K a) (4) 
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where 
u, v, and w are the components of the velocity vector 
K is the eddy viscosity coefficient 


Studies have shown that the following approximation is possible: 
£1k oQ)| x 20 (5) 
VAL ON Oz 


where; 
the condition for neglecting the horizontal terms is that they 
approximately cancel each other. Since those terms have been 
shown to be of the same order of magnitude and opposite signs, 
this condition is possible.'' (Overstreet and Rattray 1969) 


Assuming condition (5) the canservation equation becomes: 


> oi : ah) ys al 30 
U go Se ; Sle x) “Bul a) © 


Assuming isotropy and horizontal homogeneity the two dimensional 


Equation (6) becomes: 


40 9202 p[FOQFO) 
U sx er = k,| So a0 | (7) 


or 


ee ol), P20 - 20 2 2Q (8) 
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where Pp = - 
Lb 


P can be computed for each grid point if the direction of the current 
at each grid point is known. However, unless direct measurements of 
the current have been made near each grid point this 1s not generally 
the case. Therefore, the next best solution is to compute the current 
for each station using available data and extrapolate or interpolate 
it to the nearest grid points assuming that the current is uniform in 
the area. 


The current consists of four components 


l 
! 


V- 


Sh 
+ 
aX | 
> 
= 


(10) 


V, is the drift current due to wind action 
Vs is~the geostrophic current computed from the density field 
Uo is the barotropic current 
VY is the tide current 
If a strong barotropic current, vi , is present in the region to be 
analyzed, climatological data would have to be considered for its 
definition. 

Tidal currents may be significant because oceanographic data 
is only quasi-synoptic as noted before. However, currents considered 
are generally time averaged for the entire time of all the observations. 


Therefore, observations would represent a time lapse of two or three 


days or at least 24 hours. It is possible, therefore, to neglect tidal 
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currents in the average. However, when the observations are truly 
synoptic, tidal currents must be considered. i. &- 


"pidal currents are simply the water movements generated by tidal 
forces. Since tidal forces which produce tidal currents are peri- 
odical and these currents are generally only important in shallow 
water, tidal currents will not be considered as part of the ad- 
vection problem." (Ocean Thermal Structure Forecasting, ASWEPS 
manual series, V- 5, p. 21) 


Hereafter we neglect tidal and barotropic currents SO the 


horizontal components of the current vector become: 


U U 4 a5 Us 


i 


(11) 


w= Vat yY 


and 
- Verh 
Uz + Uy 
Similar formulae have been suggested by Arthur (1966) and Adem 
mo70). -— 


These velocity components can be computed from the data avail- 
able in a standard oceanographic station. 

Considering the drift current, we examine possible relation 
between it and the obsereved wind. Denoting the surface wind velocity 


as Va, the following formula apply for Va along the Y-axis: 


_(%)z 


Ui= VE Cos( 45°- F 2) 


(12) 


_(%)2 
WE Aim (4s zz) 
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Va 


2) 








where 


V/ = 0.0253 }l Va 


(Ain 
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Te 3.67 4/l 
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for WAS 6 pn see. 


and 
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(p is the latitude 
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For the general case, when the Y-axis is not along the wind direction 
Equations (12) becomes (Adem 1970), 


ae 


- 0.0126 (U,CAO-Vas4m?é 
U, = gies (1, Co | 


(13) 
Wes 2.9126 (Y, Cor O - U, tum é | 
d ; a 
sin 

Where G is the angle between the surface current and the surface 
wind. o is positive when the current is to the right of the 
surface wind. 

The geostrophic current ( Ue, Ve ), between two stations 


may be computed from the dynamic height anomaly using the Ekman 


approximation: 


a 





! gh. - nad - 79 2AD 
f yey = 4 hy 10 5y VO Oy 
2 24D - 
a Di igs A 
4 Up 4 Px = 1) EX srenely Ox 


where i is the Coriolis parameter 
AD is the geopotential anomaly 
Assuming that the density increases steadily with depth, the relative 


geostrophic current between two stations A and B, can be expressed as 


Mz Ye = BEAD.” APs) = 


where L is the distance between A and B 
A further simplification of Equation (11) may be achieved 
assuming that the geostrophic current has no contribution to the 
advection of density anamolies; 
the major contribution to advection (of heat) comes from the 
wind drift current since the slop current is assumed to be 
geostrophic, or flowing rarallel to the isotherms. That is, 
there is no advection (of heat) by water flowing parallel to the 


isotherms; only cross section flow causes a change in the tem- 


perture at a point." (Ocean Thermal Structure Forecast, SP-105, 
ASWEPS manual series, v. 5, p. 33) 


If this approximation can be accepted, Equation (11) becomes, 


2 Va 


(16) 
eel 
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P= We 
Ua 

Therefore, many degrees of approximation may be used to compute 
P and they can be summarized in three cases: 

a. Currents at the station are known. If the true current has 
been measured at each station, P may be computed directly or if the 
current has been computed, P may be computed from Equation GOT, dale 
same P may be used to nearby grid points provided the changes in current 
direction are not too large relative to the grid scale. 

b. The wind velocity is known but the true current _is not 
known for each station. In this case only the drift current can be 
considered and Equation (16) must be used to obtain the direction of the 
current for the grid points using the wind direction nearby. Also, 
in this case empirical consideration may be used for allowing fer the 
drift current. These empirical considerations are: When slope and 
gradient currents (due to mass transport by the wind) have been set up; 
they affect the direction of the original wind driven current, and when 
a steady state has been reached, the deflection of the surface from 
the down wind direction will be reduced to near 20 degrees at low 
wind speeds or to less than 12 degrees at high wind speeds (Summary of 
Technical Reports of the National Research Committee, v. 6A), (Ocean 
Thermal Structure Forecasting, SP-105, ASWEPS manual series, V. 5); 
(Hubert 1964). 

c. Only the wind field is known but other factor may be 
important. This situation may arise when the coast or bottom topography 


would invalidate the assumption that the current is at a constant angle 
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to the wind direction. For this case it is possible to use the follow- 
ing relation to obtain as estimate of the direction of the current at 


each grid point. 


Ke VV (Kitaigorodskii 1961) 
and therefore, 
LVL = COSTANT = RATIO 
K 


— 


where V is the fluctuating velocity vector. 
However, sufficient surrounding data are required to apply a numerical 
solution for obtaining these variables. In order to compute Lv! 
a P value, from the wind field is assumed for each grid point as in 
b. or for the whole grid if the wind field is more or less uniform. 
Using the assumed P, values of U/K are computed for all possible erid 
points using Equation (8). The numerical sheme to perform this will 
be explained in the next section. 

De fining Vv. _ us P ar 

K K K 

and assuming that eddies present in the area may introduce a great 
range of directions in the fluctuating velocity and given that LVL 
is a constant, the maximum value of U/K should correspond to a minimum 


value of VU). ACOs ge Ung is small enough, it is possible the 


approximation 


IC 


lV! v Os maximum 
K K 


RATIO 
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This approximation can be made if it is assumed that the mean current 
is thermohaline and therefore flows parallel to isotherms so it does 


not contribution in the advection of heat.1. @- 
ar eee, 
= 2 
WVT+tV VT = VT 
Ws POE 


Therefore, only the fluctuating component is obtained when vy 5 
and Ags are computed from a field of temperature to solve for 
V . For each grid point at which U/K is computed a’U/K can be 


obtained using the formula 


1, 2 Y, 


(Heo) © 


- : 


i 


—— 


rales 


A 


and the new P = ULK_ = 


SS 


Using this U/K and the new P Equation (8) can be solved for T at next 
eid point if enough nearby grid values exist to use the rotational 
computing molecule presented in the numerical model (seccion 2). 

The procedure above can be repeated many times but it has been 
observed that usually the first approximation is satisfactory. Because 
the quadratic equation (17) has two solutions, the wind field is used 
in the selection of the correct solution. This choice can be objective 


and included in the computer program. 
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2. The Numerical Model 


Knowing P, U/K can be computed whenever the parameter is 
known in four neighboring grid points and then, using this value of 
U/K and the corresponding P the next unknown grid point parameter can 


be computed. Tne problem reduces to solving the differential equations 


uf O +P 20 _ 20+ 2 
K, Ox 04 “oy? (8) 


- 20 a 


>|> 
= 


The corresponding Finite Differential Equations may be solved using 
the same scheme. 

The following is a description of a direct solution of the 
partial differential equations by means of a rotating computational 
molecule in such a way that a solution is possible even though data are 
scarce. 

The horizontal continuity equation, Equation (8), may be 


approximated by the following computational molecules: 


(1) a. 


Z, (Kt) = 4.08 ZK, L412) -Z (KL92)~ 2 (Ke2,L)-Z (2, 42) 


ra Be 
ae (y+ } Iz cay +4)~ 2 (w02, +2) 


ei} 





flees Z, (Kt) - (ur)? [2 (xs2)- 2 (x12) 


20x KxXR 


(2a Z (44) > ~ (236 Jlzcrnsr -20c8 casa ez orsted| 


+2.0%2 (KLt1)-Z2(K4t2)- Pu ll = | 
py ra ) Pall, [Ztss+0 Z (KLt2) 





isc Zlws)2 Zt) ~ 5 [ 2-427 -Z044 a 


Ka xXSc 


(3) a. HALES ib) ee oO Z (k-1141)- 2 (k,Lt2) ~ZK2 1) -2 (k-2, 172) 


_ (uy [2 (k-4, 241) ~2 (4-2 z,1+2)| 


Kaxk i 


Z (ki) = ZK4) - (urn) [242 2 (2,4) 


Da- 2.0"*K* KR 


(4) a. = vee Lz Cey tt) 2082 LI EOCLED| 


+ 2.0*%2 (k-4L)-2(K-24 ce ; k-tL)~ e 
) rT EE Ey) 2 yl) 


: . ae Be Ye (eos 
De Z (kDa LK SE |24 La Caray 2)| 
(5) a. 4 (K4) = fn EKG bt) BKB L)- 2K LA) A(R 2) 
P , 
= (usar) [2 (ge) 222) 


Kx XR 
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ACS AA = (uss ar%) |Z (12) -Z(K2,4)| 


2.0*K» XR 
(6) a. ZH) _ zal (K+t, L-1)-2.0*2Z (Kl-1)+ Z (K-50-9) 
4 20° Z(kuct) - B(Kl-2)- Cyd zak] 
be (kt) = Z(k1)- pH —a\z (Ke41-1) - 2 crue) 
ee i (CI) es TE eh (en ep) a(S IL) = AGG IER) A Gia) 
- Liceary" |Z (K+I t-1) — Z (2, 2) | 


b. ZA (KL) = Z, (KL) — ail 2 (41-0) 


27 7 
= 


(8) ae f is) _ / aiGh (Kei Let }- 2.082 (K+i, i) + Bite, é L); 


eeiows 


= 4+ 2.042 (KitL)- 2 (k+24)— 4 | 2 (we4e)- 2 (R20) 
Kaxse 


b. £ (kz) = 2-3 E GEES | L058 1) am) 


(9) a. a (Sie __Zcne41) 188) th) 2 (Ke + OeLO| 
[a0-(H85) +2.0] 





Z (k,4) = Z (kt) ~[2 (+44) -2 (K-44) | waz au 


[z (x,4+1) ~z (k, t-1)| 


env A iliaaae 
2.0" K* YSC 


ag 





In each case the computing molecule has been divided in two parts in 
order to give more flexibility to the computation and be able to compute 
the Lapacean and the gradient separately. In order to be able to use 
all molecules, the X-axis scale, XSC, must be equal to the Y-axis 

scale, YSC. If the preceeding is not the case, local distortion may be 
introduced. Also the above procedures are not iterative and, therefore, 


stability problem will not arise. 
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FIRGURE 4. Sequence of rotation of the computing mole: ile around 
a grid point with no value U(K,L). 
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FIGURE 5. A complete rotation of the computing molecv‘e 
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III. DESCRIPTION OF DATA 

In order to verify the proposed scheme, data were selected from 
areas with good data coverage and well known physical features. 

One area considered was the California coast between San Francisco 
Bay and Monterey Bay and the corresponding data were obtained by the 
AcOR* cruises of September 1969 and 28 April - 4 May 1970. These data 
correspond to standard oceanographic stations (Figure 2). However, 
not all the information was readily available due to the fact that the 
data were reduced for specific uses. The data available were; 

a. Wind 

b. depth 

c. temperature 

d. Sound velocity 

e. Beam transmittance 
at observation depths only. Under these circumstances the general 
procedure, Equation (11), doesn't apply and therefore. the modified 
procedure, c., was used. 

A second area was considered in order to verify the objective 
analysis procedures at low latitudes. This second set of data was 
obtained by the Colombian Navy's Division of oceanography in the 


Caribbean sea off the Colombian coast. In this case all the data from 


1 
Cruises by ships of the U. S. Oceanographic Office under inves- 
tigation plans by the Naval Postgraduate School. 
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a standard oceanographic Station were available. However, the same 


simplified procedure, C., was used in order to obtain results com- 


parable to those for the California coasts case. 


e) 9) 





IV. DISCUSSION OF RESULTS 


A. GOMPARISON OF OBJECTIVE AND SUBJECTIVE ANALYSIS 

A subjective analysis, Shepard (1970), of sea surface temperature 
for the area along the California coast appears in Figure 6. Shepard 
performed the analysis in order to make a comparison of oceanographic 
parameters during upwelling. The objective analysis of the same data 
appears in Figure Ts 

For the objective analysis, Figure 7, the grid was rotated 47 degrees 
in order to follow the general trend of the coast. Because the mean 
wind was generally along the X-axis during the cruise, an invita 
value of P equal -0.27 or a current 20 degrees to the right of the 
X-axis was selected. 

A comparison clearly shows greater complexity and more details in 
the objective analysis than the subjective analysis. The obvious 
question that arises is; do these details represent artificial additions 
or real features with physical meaning? Also, do they represent features 
expected with the inherent assumption of the objective scheme? 

It may be stated that subjective analysis is a linear interpolation 
which normally depends only on the available data for one oceanographic 
variable at a time. This generally excludes consideration of currents. 
Representations in intermediate regions with no data are generally 
smoothed contours. The resulting scalar field (as depicted in Figure 6) 
is very familiar in the oceanographic work and a departure from this 


picture is often viewed with suspicion. Adverse appraisal of a more 
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detailed representation is not totally justified as stated in the 
following passage: 
"Linearity between widely scattered data should not be assumed. 
It is a paradox of analysis that the less data the analyst has, 
the more the analysis show smooth curves and evenly spaced per- 
turbations. It is only when more data are added and the presence 
of detailed structures becomes apparent but not entirely delineated 
by the data that the analysis becomes difficult.'’ (Ocean Thermal 
Structure Forecasting, ASWEPS anual Gerveceevee >, Ds 27m 

Other features present in the objective and absent in the subjective 
analysis are the circular patterns or pools of relative warmer or cooler 
water. This type of temperature distribution is intuitively locical 
as most swimmers can tell and physically may be expected. Subjective 
analysis procedures are not suitable for detecting them because they 
do not normally involve considerations of temperature along with 
advection currents and diffusion terms. 

''rddies have seldom been fully delineated. Usually a single line 
of temperature observations crosses what may be a roughly circular 
core of warmer water or a tongue of a warmer water extending back 
to the source." (Summary of Technical Reports of the National 
Defense Research Committee, V- 6A, p. 73) 
The importance of including the influence of the wind action, i. é. 
wind currents, within the scalar analysis is emphasized further in the 
following physical arguments. 

Wind and wave action have a net result of mass transport of water 
and therefore both are capable of producing zones of convergence or 
divergence accompanied by a rising or lowering of the thermocline or 
may be responsible for water of one temperature being carried over 
water of another. (Summary of Technical Reprots of the National 
Defense Research Committee, vV- 6A). For example, 

the effect of the wind in maintaining or deepening an isothermal 


surface layer is complicated by the fact that the wind setsS up 4 
current near the surface which may remove the mixed .ater from a 
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given area almost as fase AS 2b Bomeommea— = (Summary of Technical 
Reports of the National Defense Research Committee, V- 6A, p. 42) 
B. INTERPRETATION OF RESULTS IN TERMS OF THE ANALYTICAL MODEL. 

Additional interpretations of the features within the objective 
analysis are possible by considering the basic processes included in 
the analytical model. The physical description of the analytical model, 
Equation (4), used in the objective analysis scheme is that when a steady 
state is reached the advection of the scalar property is equal to the 
eddy diffusion of the scalar in the medium. With this physical inter- 
pretation, we can explain certain features of the circular distribution 
of properties when these features are frozen in time (stationary). 

These features are: 

1. Circular patterns, 

2. Tongue type patterns, 

3. Gradients 

4. Size of closed isolines. 

With respect to the circular patterns, heating and cooling of the 
sea surface is not uniform for the entire layer of the analysis and 
therefore in any event a local circulation will be generated in order 
to obtain the equilibrium of mass. Warm water will flow from the 
source toward the periphery in all directions and the Coriolis accel- 
eration will deflect the paths to the right (in the Northern Hemisphere ) 
until they become circular. This situation will exist as long as a 
disequilibrium in mass exists and will become a steady state when 
advection out of the circle is equal to the inflow of warm water, 4 
similar situation is reached in a cold seal where warm water from the 


periphery is going to the center of the pool following a circular path. 


41 





Decaying drift currents will also produce the same effect because once 
the wind ceases the drift current continues to move for a certain time 
under its inertia while being acted upon by the Coriolis acceleration 
and while being damped by the diffussion term. 

With respect to tongue type features, the oceanographic literature 
has abundant examples of tongue-like diffusion of properties such as. 
salinity and temperature. For example, it has been observed that in 
areas of convergence, convection cells, the surface water moves down- 
wind at a considerably more rapid rate than in areas of divergence and 
therefore the flow of water in response to the wind should be uneven. 

With respect to the gradients, because vertical mass transport 
occurs in the convective layer, it is clear that the assumption of a 


horizontally homogeneous mixed layer above the thermocline is an over- 


oR 


simplified concept and relativeiy high graaients may be possibile in 


mf 


{p 
eo) 


short horizontal distances. This is especially true in the area 
examined which is in coastal waters, where an irregular wind field 

may have existed and where upwelling is known to exist. The analytical 
model was, of course, a non divergent and this may not have been the 
case. Therefore, greater than normal gradients inside the inertial 
circles may reflect over compensation for the upwelling or downwe Lling 
effect in the temperature. Although, this is, perhaps, an unrealistic 
result, it could be smoothed out by computing the amount of this effect 
using the grid points over nearby coastal regions (imo tect Numerical 
Weather Facility. Technical memo No. 5, January 1965). On the othe 
hand, artificial intensification of the gradient may still be used to 
identify convergence and divergence Zones by allowing a maximum gradient 


in a given region depending of local conditions. 
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With respect to the size of closed isolines, the radius of the 
inertial circle for the type of circulation Eaaiaioed above varies 
with the velocity and the sine of the latitude: For a current of 
average velocity located in mid-latitudes this radius is about 15 miles 
which can be greater if the Pereulation 1s witimenced by va horizontal 
component of gravitational force or other lateral force (Summary of 
Technical Reprots of the National Defense Research Committee, v- 6A). 
It is interesting to note that 15 miles is the order of magnitude 
observed in the circular features of the objective analysis under 
consideration. 

Therefore, it seems possible, with this analytical model, to relate 
circular and elliptical features to jnertial currents due to differ- 
ential heating and cooling, convergences, divergences, decaying wind 
drift currents; tongue like features to normal diffusion processes, 
and intensified gradients inside circular patterns to strong upwelling 
or downwelfing. However, gradients are represented very well within 


the analytical model. Reasons for the latter feature are however, known. 


C. COMPARISON OF OBJECTIVE ANALYSIS RESULTS FOR DIFFERENT LEVELS 

In order to verify vertical consistency, a mixed layer of at tease 
10 meters was assumed and an objective analysis of the 10-meter level 
was made and is shown in Figure 9. However, the wind field was not 
uniform and, there was not enongh information about fetch and duration 
with available wind observations. Therefore, we do not expect to obtain 
an exact correspondence of the features of the two layers over Ene 
entire grid. Nevertheless, vertical continuity is observed when the 


main features of the two analyses are superimposed in “igure 10. 
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Maximum and minimum temperature regions are marked with letters A 
through K in Figure 7 and apparent corresponding features in Figure 9 
have been marked with same letters with a prime sign. The same regions 
are shown in Figure 10 but delineated by only one or two isotherms in 
order to simplify the resulting picture. 

The comparison of thermal structure at the two levels shows an 
almost identical correspondence between regions 

A and A‘ 

B and B' 

D and D' 

E and E' 

F and F' 

H and H' 

T and i' 

J anda. 
This may be expected from good vertical mixing and more or less constant 
advection in the entire layer. 

Regions CeandeG: 

K and K' 

G and G' 
have horizontal relative displacement from 5 to 12 miles. This may be 
due to a different velocity field for the two levels. Past history 
(i. e. two or three weeks earlier) would show how these advective 
distortions were generated. Similar results may be obtained for the 
September data (Figures jiltmand.= 12 )- 

Considering that each analysis is completely independent from the 
other, this comparison shows convincing verification ef Verelucal 


consistence. 
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D. OBJECTIVE ANALYSES OF SURFACE TEMPERATURE, SALINITY AND SIGMA-T 

FOR LOW LATITUDES 

Objective analyses of temperature, salinity and sigma-t for a region 
with mean latitude of 10° north appear in Figures 14, 15) fande 67 ie 
main reason for considering these analyses for a region not very well 
studied from the oceanographic point of view is to compare the effect 
of latitude in the size of those features interpreted to be inertial 
circles in the preceding discussions. 

The region of California coast considered has a mean latitude of 
37° north and feature A in Figure 7 appears to be representative one fOr 
this region. It is observed that its mean radius is about 8 nautical 
miles. The radius of inertial circles is inversely proportional to the 
sine of the latitude and therefore it would be expected that the radius 
of the inertial circle be increased by a factor of approximately 3.4 
when translated from latitude 37°N. to 10°N. a corresponding feature 
in Figure 14 was marked A' has mean radius of about 35 nautical miles. 
Therefore, the test for latitude effect seems satisfactory from this 
gross comparison. However, more precise computations would be required 
before a conclusive result may be obtained. 

The region off the Colombian coast, considered in the low latitude 
study (Figure 13), receives the run-off of three main rivers: ® the 
Atrato in the south-west region a the grid and the Sinu and Magdalena 
rivers in the east. Features L,M and N in Figure 15 correspond to places 
where the run off occurs. The patterns observed at these points coin- 
cide very well with what may be expected and are very well correlated 


with the sigma-t distribution, L', M' and N' in Figure 16. 
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V. CONCLUSIONS AND SUGGESTION FOR FUTURE RESEARCH 


A. CONCLUSIONS 

A small scale, high resolution and low computing time objective 
analysis model has been developed. More tests and field verification 
are required to show all of its advantages and shortcomings but tests 
described and interpreted in this study suggest that if has great 
possibilities to become a basic step in forecasting oceanographic 
parameters on a tactical unit level. One desirable feature OF eae 
results was the detail obtained in contrast with that obtained in a 
subjective analysis. The required input and simplicity enables it to 


be used to obtain an almost inStantaneous analysis during an oceanographic 


Fh 


survey and therefore to investigate interesting features mot soreseen 
in planning stages. Therefore, it could be applied for the purposes of 
reducing costs and increasing scientific results of an oceanographic 
expedition. 

The procedures described need to be refined in order to allow for 
convergence or divergence whenever these phenomena become important. 
However, by understanding how the analytical model (section B) handles 
these features, meaningful interpretations can be made even for these 
features in the results. The data should be corrected for diurnal 
effects but there is not a satisfactory procedure for this (Technical 
report, ASWEPS No. 13). An alternative solution, observations at the 


same hour of the day, may be undertaken by XBT drops. This is perhaps, 


unfeasible at most oceanographic stations due to operational costs. 
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B. SUGGESTIONS FOR FUTURE RESEARCH 

A suggestion for future investigation would be to integrate this 
scheme into a forecast procedure. The following steps could be applied 
using the present form of the scheme: 

a. Objective analysis using horizontal grid at standard depths. 

b. Conversion to a three-dimensional grid using the vertical grid 
as a complementary tool. 

c. Repetition of procedures 4. and b. at equally spaced time steps. 

d. Determination of Lagrangian time rate of change for identificable. 
structures (assuming that some of them are persistent in relation to the 
time step selected as seems to be the case in the neighborhood of 
constant currents). 


e. Determination of changes in the heat budget for each character- 


whi 


istic structure and ¢elating) it to easily measured parameters Seon 
amount of clouds in each Standard level, sun's declination (season) and 
hour of the day. 

£. Measurement of advection by means of changes in.the structures 
in point d. and (relating) it to permanent currents, wind speed, direction 
and duration. 

g. Forecast for the next time step. 

h. Forecast verification and obtention of corrective cefficients. 

In the delineation af the steps leading to a forecast it has been 
assumed that some features may have a relative long duration or per™ 
sistence but this situation must be investigated for each application. 

If it is proved, identification of persistent or semipersistent features, 
their changes in time and parameters affecting them may be a simp le 
approach for detailed forecasts, even for handling treasient thermocline» 


now practically impossible to forecast. 
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